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can be used to measure the spin of the invisible graviton, or of any other invisible objects 
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confirm, especially for the case of two extra dimensions, that the method is valid within 
a realistic setup, which includes the simulation of the Standard Model background, beam- 
strahlung, bremsstrahlung, calorimeter resolution and calorimeter coverage. When the 
number of extra dimensions is increased, the angular distribution does not provide suffi- 
cient information to extract the spin, but this method still offers a useful parameterization 
of the single photon cross section using which the nature of the missing object can be 
studied. 
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1. Introduction 

1.1 Theory Background 

The Standard Model (SM) of Particle Physics, which has so far stood the test of numerous 
precision measurements, is nevertheless not considered to be the ultimate theory. One of 
the reasons is a serious shortcoming of the SM, namely that the mass of the Higgs boson 
diverges quadratically as we evolve the energy scale of the theory towards the electroweak 
scale from the Planck scale. This is the so-called Hierarchy Problem. 

We expect that the SM as an effective theory is valid up to 0(1 TeV), and any signs 
of New Physics, if they should exist, will arise at that energy scale. Supersymmetry is 
a well-motivated theory which solves the Hierarchy Problem, but another possibility has 
been proposed [|| recently. According to this theory, the SM particles are confined in a 3- 
brane, and only gravity propagates in the compactified extra dimensions. The weakness of 
gravity is explained by the fact that it propagates in the extra dimensions, and the overlap 
of the wavefunctions between the graviton and the SM particles is small. The fundamental 
scale of the theory is not Mpi ~ 2.4 x 10^^ GeV, but Md = 0(1 TeV). Here: 

m2 ~ KM^+^ (1.1) 

where n is the number of extra dimensions and Vn is the volume of the extra dimensions. 
Thus the weak scale and the fundamental scale are very close, and so the Hierarchy 
Problem is 'resolved'. It may be argued that this 'resolution' only translates one problem, 
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namely the hierarchy between the electroweak scale and the fundamental scale, into an- 
other, namely the size of the extra dimensions. Let us nevertheless adopt this model as an 
effective theory applicable at low energies, leading to concrete predictions which should be 
studied. 

Since the graviton propagates in the compactified extra dimensions, there are many 
massive Kaluza-Klein excitations. Such excitations couple to the SM particles, and we can 
expect to observe at future colliders massive graviton excitations, which we call massive 
gravitons hereafter for simplicity. 

Some bounds are available already from current collider experiments. DO Q and L3 
Q have made searches both in the direct production channel where the produced graviton 
is undetected and is identified through the missing 4-momenta, and in the indirect channel 
where virtual graviton contributions give rise to 'contact term' type effective interactions. 
The direct limit on M^, can be set only in the former of these two search modes, and the 
current limit from L3 |^] is Md ^ 1 TeV for two extra dimensions. 

When the number of extra dimensions is low, especially when n = 2, Astrophysics sets 
a stringent bound on the fundamental scale M/j. Graviton emission into extra dimensions 
from a hot supernova core is considered in Ref. |^], and the phenomenology of SN1987A 
places a strong constraint on this energy loss mechanism, allowing us to derive a bound on 
the fundamental scale M^). When n = 2, the allowed lower bound on the fundamental scale 
M/) is Md > 50 TeV. Furthermore, the graviton decay contribution to the cosmic diffuse 
gamma ray radiation sets the bound M^, ^110 TeV for the case n = 2 j^. The bound 
becomes much less stringent if we raise the number of extra dimensions n. The constraint 
from SN1987A becomes Mq > 4 TeV for n = 3, and Mo > 1 TeV for n = 4. There is 
theoretical uncertainty associated with these numbers, and so the TeV fundamental scale 
is not ruled out. In any case it is meaningful to derive alternative limits and so we do not 
limit ourselves to parameter regions not excluded by Astrophysics. 

General studies of graviton measurement have been made in 0, § • Search strategies 
for future e^e~ colliders the LHC [|ll|, the 77 colliders [Q, and 67 colliders |jl3[ 

have also been studied. Recently, there has been increased attention on the higher energy 
blackhole and 'trans-Planckian' physics which may be observed at LHC [0]. We do not 
consider these effects as the energy scale is lower in our analysis, and hence the leading extra 
dimensional effects, in the form of Kaluza-Klein excitations of the graviton, are expected 
to be perturbatively calculable. 

1.2 Motivation 

We note that the above studies have concentrated mainly on the discovery of the extra 
dimensional phenomenology. In this regard, although the difference between the signal 
and the SM background has certainly been considered and exploited in the above studies 



|6[-|13], no systematic study has been available so far, for the case of these 'invisible' 
gravitons, to claim the discovery of spin-2 massive excitations. 

Let us discuss this point further. For the case of massive spin-2 resonances, the spin-2 
nature is probed by measuring the decay angle 0. The presence of terms proportional to 
cos^O then provides evidence for the spin-2 nature of the excitation. 
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This same measurement can be made with the angular distribution resulting from vir- 
tual graviton exchange. One may measure processes such as e~e^ — > 77 and the presence 
of terms proportional to cos^0 can be confirmed in principle. This, however, is not con- 
clusive evidence of spin-2 gravitons, as all that one is measuring here is the nature of the 
'contact-term' type interactions, L cx Tsm'^'^Tsm^i/, which arise from integrating out virtual 
graviton contributions. 

In order to claim the discovery of spin-2 massive excitations it is necessary, therefore, 
to produce them directly, and measure distributions associated with them. This is the 
topic of our present exposition. 

As each of the massive graviton modes has a coupling suppressed by M^;, it escapes 
detection and the events are characterized by missing 4- momenta. At hadron colliders, only 
the missing transverse momentum is measurable, and so the extraction of the spin is, even if 
this is possible at all, expected to be highly model dependent. Let us therefore concentrate 
on the case of lepton collisions, more specifically the e~ future linear colliders. 

At first sight, it may appear that the measurement of polarizations participating in the 
interaction is sufficient. This is unfortunately not so. For e^e"*" annihilation the polarization 
structure is exactly the same for both spin-1 and spin-2 objects. For 77 annihilation the 
anticorrelation of the photon helicities is an interesting property of the direct production 
of spin-2 objects, but it is not possible to trigger such an event as no measurable particle 
is produced. Even in 67 collisions it is not possible to utilize only the polarizations to 
determine the spin of the excitation. 

This leads us to the consideration of the process e~ G7 as the only way to 
measure the 'spin' of the graviton excitation. 

There are two variables which we can measure. These are the energy and the polar 
angle of the photon. The energy of the photon measures the mass of the excitation. A 
study is available in [|lO| which compares the distribution of graviton events against that 
of the supersymmetric e~ XiXiT- Our approach is more systematic and rigourous. 

Our starting point is the realization that the collinear divergence, which dominates 
the G7 cross section, is only logarithmic. At the amplitude squared level, the divergence 
scales as l/p^^ where pr, or Et which is the quantity that is actually measured, is the 
photon transverse momentum. Thus it is possible to regularize the collinear divergence of 
any 7-l-missing events by multiplying the distribution by or Elj.. 

The resulting distribution, which we may call the E'^ moment distribution, shows a 
remarkable behaviour. In short, spin-5 excitations have moment distributions which con- 
tain terms up to cos^'^0, or equivalently cos(2S'0). Hence by Fourier analysis, in principle, 
we can extract the nature of the excitation. 

The SM background process which is at the leading order in the electroweak coupling 
constant, e~ — > Z^^ — > vv^^ only contains terms up to cos 2^. Furthermore, we can 
reduce this background by cuts around the corresponding photon energy. 

The case of the subleading order background process, e~e^ — > I'e'^el, which occurs both 
via virtual exchange and via the off-shell Z^, is less straightforward. We also calculate 
this background. The induced contribution can contain terms proportional to cos 40 
and higher. Hence this background must be subtracted from the measured distribution. 
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As the analysis relies on measuring only one visible object, it is necessary to carefully 
consider the systematics which can affect the measurement. In addition to the calorimeter 
resolution and the angular coverage, there are further obstacles from initial state radiation 
(ISR), composed of bremsstrahlung and beamstrahlung. 

We concentrate on analysing the systematics in this paper, and leave the statistical 
analysis and the error estimation, as well as a more detailed account of the hard/multiple 
photon emission effect, as topics for future studies. In this regard, we should note that the 
analysis presented in this paper is realistic, but the numbers are not conclusive. 

We finally note that the work presented in this paper provides a general framework for 
analysing events composed of 7+missing 4-momenta, which can be used to search for any 
New Physics, not only those associated with extra spacial dimensions. As an example we 
discuss the case of the supersymmetric XiXil process. 

This paper is organized as follows: in Sect. ^ we discuss our approach and explain the 
procedures adopted to make realistic simulation. In Sect. ^ we present the results of our 
calculation. The conclusions are stated at the end in Sect. |^. 



2. Method and Discussions 



2.1 Fourier Analysis 

Let us first define the kinematic variables 6 and x^. 9 is the angle between the photon 
and the momentum direction of the incoming electron. There is no forward-backward 
asymmetry and so our results are unchanged if we measure the angle from the positron 
direction. The energy fraction is defined by: 

^- ,2.1) 



y/ '^nominal 



In the following, s = Snominai denotes the nominal centre-of-mass energy squared. The 
effective centre-of-mass energy squared is written as s. The doubly differential distribution, 
(P a / d cos 9dx^, has both soft and collinear divergences, and therefore we expect that there 
is little visible distinction between the signal distribution and the background distribution. 
For the sake of clarification, for the two-body case, disregarding the effect of ISR, we have: 



d cos 9dx-y 327rs 



|A^|2, (2.2) 



where is the corresponding matrix element squared summed over final state helicities 
and averaged over initial state helicities, and for the three-body case we have: 

d'^a x^ f di^ 



■\M\^. (2.3) 



dcos9dx^ 512-11^ J Att 

The integration in the three-body case is the integration over the solid angle of the invisible 
'decay' in the rest frame of the 'decaying' system. 
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As stated in the introduction, this divergence is only logarithmic, and so we can elimi- 
nate this divergence by multiplying the distribution by the transverse energy squared, E^, 
of the photon, viz: 

£ < a El > EldPa 

^ — ^ (2.4) 



dcosOdx^ d cos 6dx^ 
We note that an alternative definition, which differs from the above definition by factor 
E^, may be more useful in experimental analyses: 

d'^ < gEt > _ Erd^a 
dOdx^ dOdx^ 

We take the first definition as we consider it theoretically more elegant. 

This E"^ moment distribution has the empirical property, which we have verified an- 
alytically for several cases and we attempt to justify in the next Subsection, that only 
angular dependencies up to cos^^O are present for the production of invisible spin-5 ob- 
jects. Assuming this property, and noting that cos°^9 can be expressed in terms of a sum 
of terms proportional to cos/30, (3 < a, we can apply the Fourier expansion to this moment 
distribution. 



E^d^a 
d cos 9dx^ 



ao(x^) + ai{x^) cos 29 + a2{x^) cos 49 . . . 



s 



= ""^^ amix-y) cos 2m9 . (2-6) 

m=0 

The coefficients obtained, at least in principle, as: 

am{x-y) = - —^——cos2m9d{29) (m = 1, 2, 3, . . .). (2.8) 
vr Jo d cos 9dxy 

Hence the measurement of the maximum m for which is nonvanishing gives a measure 
of the 'spin' of the produced object. We note that the dimensionless quantities, 

which we measure in units of fbGeV^. Fourier analysis, as presented here, is not the only 
way to extract these coefficients, but it is the most theoretically elegant. It can also be 
considered, at the very least, as a way of parameterizing the two dimensional distribution 
of Xy and cos 9. 

In the presence of ISR, when evaluating the Fourier coefficients, we need to be more 
careful as to which frame x^ and cos 9 are measured in. To be consistent, these should be 
evaluated in the lab frame throughout. However, it is more convenient to carry out the 
calculation in the centre-of-mass frame using the Monte-Carlo approach. In this case, one 
must be careful to include the correct Jacobian factors. 

We have tabulated the leading order theoretical values for these coefficients in Tabs. || 
and The coefficients are obtained by multiplying the functions in Tab. § by the prefactors 
in Tab. |l] and h"^ = 0.389 x 10^^ fbGeV^. Further discussions are given in the Subsections 
to follow. 
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Table 1: Overall prefactor functions for the Fourier expansion. One should also multiply by 
= 0.389 X 10^^ fl) GcV^. z^ei^e in the above table refers to the exchange, in the limit of zero 
momentum transfer along the virtual propagator, 'spin 1' and 'spin 0' tabulated above arc for 
the production of heavy bosons followed by the decay to two massless particles. The photino pair 
production corresponds to the limit of heavy selectrons. 
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Table 2: Coefficient functions for the Fourier expansion. 



2.2 Spin and Angular Distribution 

The naive technical reason for the correlation between the spin S and the number of 

coefficients mmax of the Fourier expansion is as follows. 

A spin S boson is a rank S tensor with S indices. Each of these indices, when contracted 
with the electron current or a momentum operator, gives a contribution to the amplitude 
which is proportional to up to cos0. Hence we get cos^'^^ at the amplitude squared level. 

Perhaps a more physical restatement is that the plane wave due to the produced spin-5, 
or rank S tensorial, object intersects the electron current S times at angle 9. 

There is a loophole in this argument concerning the case of the s-channel photon 
mediated contributions. In this case the intermediate photon, which is a vector, carries the 
angular information and so the distribution is proportional to up to cos^^ already before 
multiplying by E^. Hence the moment distribution {E^cf a/ d cos 9dx^) contains terms up 
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to cos^^, regardless of the spin of the produced object. 

This is a Umitation of our method, but there are only few cases where this becomes 
problematic. If the s-channel graph is solely responsible for the production, the collinear 
singularity is absent after background subtraction, and so it becomes immediately clear that 
there is no f-channel contribution. The ambiguity concerning the 'spin' of the produced 
object arises only in cases where the s-channel and i-channel terms co-exist, and only for 
spin lower than 2, i.e., spin-0 or spin-1. 

Out of these two cases, the spin-1 case is genuinely problematic and the signal is not 
distinguishable from a spin-2 signal. However, we note that a photon-photon-vector vertex 
lacks theoretical motivation. 

If the produced object is spin-0, the s-channel term is distinguishable from the t-channel 
term by beam polarization. Hence this case is still not serious. 

Having explained this, it is instructive to now look more closely at spin-2 production, 
as this case involves both the s-channel, t-channel and contact terms. 

Out of these, the s-channel term is separately gauge invariant with respect to the 
photon, and so is the sum of the t-channcl and contact terms. We need to consider whether 
the cos*^^ behaviour originates in the latter of these contributions. Indeed, in the Feynman 
gauge, raising the tensor rank of the field only amounts to further contractions with the 
metric tensor g^^,, and so by dimension counting, and by forward-backward symmetry, we 
see that the contribution in the spin-2 case only goes up to cos^ 0. 

On the other hand, the longitudinal contributions from terms that are proportional to 
PfiPu do have the right dimensions for a cos^0 behaviour, but we generally expect that such 
contributions, coming from Goldstone modes, behave like objects with lower spin. 

Explicit evaluation shows that such contributions do lead to cos*^^ behaviour. How- 
ever, this behaviour comes from an unphysical source, namely the non-conservation of the 
energy-momentum tensor. The inclusion of the s-channel graph is necessary. This cancels 
the unphysical cos^O behaviour and leads to the final distribution which is physical and 
contains terms up to cos^0. 

It is possible to restate the preceding discussion in terms of the individual helicity 
components. However, in our opinion, this type of analysis only provides another viewpoint 
and does not enrich our insight into the nature of the problem at hand. Although it is 
clear that we have not provided a rigourous justification, to some extent, we believe that 
this is a price that we have to pay for retaining generality. 

From the above discussion, it is clear that the presence of a tensorial interaction is 
not sufficient to guarantee a cos^0 behaviour. Indeed, we have verified this for the case of 
an anomalous coupling vector boson V , whose coupling to the electron is proportional to 



and so this interaction only gives rise to a cos^^ behaviour of the E'^ moment distribution. 



{l/K)a^yFi'''. We have: 




(2.9) 
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What of the spin-3 and higher spin objects? Without exphcit example it is difficult 
to say, but we expect the above argument to hold, such that the cos^9 term can originate 
solely from an unphysical origin. 

The case where the produced object is a composite system is discussed next. 

2.3 Background and SUSY Signal 

The leading Standard Model background comes from i^v"y, which includes the contribution 
— > ui?, as well as the t-channel 1/F^-exchange leading to the final state VeT^el- The 
Feynman diagrams are as shown in Fig. |^. 




Figure 1: The Feynman diagrams for the e 
contribute to the VfiV^^i and VtVt^ final states. 



vv^ process. Only the bottom two diagrams 



From the preceding discussions, and as is easily verified by explicit calculation, the 
mediated contribution to the moment distribution is proportional to up to cos^^ only. 
It is not so obvious that this is still the case for the mediated contribution and the 
interference between the two contributions. 

The analytical expression for this process, which includes both the Z^ and con- 



tributions, is given in Eqns. (2)-(5) of Ref. [15|. There is an error in equation (3) available 
therein, and one should correct this error. The s' in the overall factor should only multiply 
the first two terms inside the square brackets. The expression without this correction is 
dimensionally inconsistent. We have performed four more independent calculations, out of 



which two are analytic and the other two use HELAS |16|, and verified that the numbers 
are consistent with each other. 

To motivate the discussion, let us first consider the mediated background in the 
limit of large mass. In this case the fusion graph (f) of Fig. |^ vanishes and we 
have only graphs (2) and (3), both of which can be written in terms of an effective contact 
interaction between two charged currents. 
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However, by Fierz transformation, this contact interaction between two charged cur- 
rents is equivalent to the contact interaction between two neutral currents, and so the 
amplitude has identical form to the mediated background [^]. Hence the distribution 
is only proportional to up to cos^^. 

Furthermore, also by Fierz transformation, the SUSY signal 777, or more generally 
XiXi7; has the same structure in the limit of large selectron mass. We thus obtain the 
expressions listed in Tabs. ||, ^. 

It is found that at LEP energies, although the above consideration of the heavy 

limit is clearly inappropriate, the shapes of angular distributions corresponding to the 

and mediated contributions are fairly similar. There is no guarantee that this 
will still hold at future LC energies, and explicit evaluation indeed shows that there is 
nonvanishing cos 40 and higher contributions. 

Even so, what becomes clear is the following. Let us consider the general case of 
the production of a photon and two invisible objects. If the particle exchanged in the t 
channel is heavy, as is the case for the SUSY XiXi7 production mediated by t-channel 
selectron exchange, the interaction is approximated by a contact interaction and then the 
distribution can only show a cos IQ behaviour. On the other hand, if the exchanged particle 
is light, as is the case for the contribution to the SM background at LC energies, the 
interaction is expected to be well understood and so the contribution can be subtracted, 
either from the bare distribution or from the Fourier coefficients. This provides partial 
justification for our procedure. 

We also note that the mediated background involves left-handed electrons and 
right-handed positrons, and therefore it can be reduced by means of beam polarization. 

2.4 Experimental Limitations 

The experimental limitations which we take into account are the calorimeter resolution, 
angular coverage and beamstrahlung. The angular resolution is not expected to be a serious 
constraint. Beamstrahlung is treated using CIRCE as discussed in the next Subsection. 
For the EM calorimeter resolution, the EM track energy resolution AE' is expected to 



be |20|, 



ono/ 

" 0.6% (TESLA). (2.11) 



V^(GeV) 

This smears the photon energy distribution. As an example, for a 100 GeV photon, Ai? 
at JLC is 1.8 GeV, wheres at TESLA AE' = 2.1 GeV. For photon energies lower than 
about 250 GeV, the energy resolution is always better at JLC. As a pessimistic estimate, 
although this choice does not affect our results in any visible way, let us adopt the TESLA 
numbers in our analysis. 

Next, let us consider the coverage of the EM calorimeter. Typical numbers for this 
coverage ||20| , ^ are in the range 50 mrad~200 mrad. In this study, let us adopt | cos0| < 
0.995, which corresponds to an angular cut-off of 6*111111 = 100 mrad pO[|. 
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The derivation of the Fourier coefficients is completely successful only when the range 
of integration is propotional to the cyclic bound. The angular coverage is expected to 
become a problem if the Fourier component under consideration changes phase near this 
unmeasured region. For a Fourier component proportional to cos(2n^), the first change 
of phase occurs at 6 = n/An. If we wish to measure n < 3, we need a coverage down to 
6*111111 < 7r/12 = 260 mrad. 

Explicit evaluation shows that the vanishing of the term a^, cos (60) is indeed contami- 
nated with the above coverage, | cos9\ < cos 6*11,111 = 0.995. In order to evade this problem, 
as a simple remedy, we perform the extrapolation, that is, the approximation that the 
values in the non-coverage regions have constant values equal to the end-points. This is 
illustrated in Fig. 0. 




original 




truncated 




extrapolated 



Figure 2: Illustrating the extrapolation procedure. The vertical axis is the E"^ moment distribu- 
tion, and the horizontal axis is 9. The dependence is not shown. 



The coefficient functions in Eqn. (p.TD and (|2.8|) are then redefined by 



00(2;^) = — 




a cos dax^ 



d{26) + 



2lT 



2[7r- 



d cos Odx^, 



(2.12) 



d{2e)\, 



E^d'^a 




cos{2me)d{2e) 
cos{2me)d{2e) 

cos{2m0)d{20) 



(2.13) 



(m= 1,2,3, 



We have found that this remedy is sufficient to overcome the problem at hand. We note 
that in the presence of ISR, the angles quoted above are measured in the laboratory frame 
rather than the centre-of-mass frame. This makes the imposition of the above cut-and- 
extrapolation procedure more difficult, and we adopt an approximate, 'wrong', solution 
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which is to impose the angular cut and extrapolation in the centre-of-mass frame. Of 
course the evaluation of the integral is otherwise performed in the laboratory frame, as 
required. 

For the centre-of-mass energy we mainly take ^/s = 500 GeV. We also consider the ^/s 
dependence by varying this to 1 TeV. We omit the discussion of statistics in this paper, 



but the luminosity may be considered to be in the region 100-500 fb per year [20, 21]. 



2.5 Initial State Radiation 

The non-singular part of the higher order perturbative contributions is expected to mainly 
affect the normalization and not the shape of the distributions. However, the logarith- 
mically enhanced bremsstrahlung corrections must be considered for realistic simulation. 
There is another source of ISR, that is, beamstrahlung. We consider both effects. 

Beamstrahlung is the radiation of photon caused by the interaction between the in- 
coming electron and positron bunches. This smears both the centre-of-mass energy and 



the momentum along the beam axis. We use CIRCE ]19| in order to estimate this effect. 

Bremsstrahlung is the radiation caused by the interaction between the electron and 
the positron participating in the annihilation event. We estimate this effect by using the 
integrated expression of |22|. This does not take into account the effect of large angle 
emission into the detector, nor does it take into account the emission of more than one 
photon from each incoming particle, but is expected to model the dominant effect well. 

As an illustration, the typical energy carried by a bremsstrahlung photon is 

/ d{E,x) (Ex) dx, (2.14) 
Jo 

where, from the expression of p^], 

d{E,x) ~ ^^(log^ - i) ~ ^ (i? ~ 250GeV). (2.15) 

TT X nip 2 x 



Hence Eqn. (|t|) yields 16 GeV for a 250 GeV beam. This energy taken away by the 



photon affects both the centre-of-mass energy and the frame. This change of frame implies 
a Lorentz boost to the final state photon and so the angular distribution is also affected, 
as well as the energy distribution. 

The subdominant 'hard bremsstrahlung', where photons enter into the detectable re- 
gion, is not treated in this paper, but we expect that a minimum Et cut is necessary to 
reduce such effect. Such a cut would also be effective in reducing the background coming 
from low Et Bhabha and forward photon events. This would be equivalent to assigning an 
energy dependent cut for the minimum polar angle ^mim and we expect that this can be 
dealt with by the extrapolation procedure adopted for calorimeter coverage in the preced- 
ing Subsection. As an example, for a 50 GeV photon, requiring Et > 10 GeV corresponds 
to 

^min — 201 mrad. Further study is desirable to quantify this effect. 
As an illustrative example, we consider a configuration in which a photon is produced 
at x^cM = l,cos0cM = 0. The Lorentz boost is easily shown to be given in terms of the 
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electron and positron momentum fractions xi and X2 by: 

Xi - X2 



cos 6liab 



Xi + X2 



(2.16) 
(2.17) 



The effective centre-of-mass energy is given by V§ = ^sxiX2- One can immediately verify 
that the transverse energy is unaffected by this boost and is given by Et = \/I/2. We 
can now plot these functions for a given centre-of-mass energy ^/s and beamstrahlung 
spectrum. Taking = 500 GeV and the TESLA beamstrahlung spectrum, we obtain the 
results shown in Fig. ^. Prom these figures, we can obtain a rough estimate of the effect of 



0.82 0.B4 0.81 



0.9 0.92 0.91 0.90 0.91 



Figure 3: Distributions of cos^? (left) and (right) of the a;^cM = 1,cos0cm — events after 
beamstrahlung and bremsstrahlung. The distributions are normalized to unity. 

ISR. cos9, or equivalently 9, is smeared typically by about 0.01 (radians). Although this 
smearing is small enough for us to expect that ISR will not change the angular distribution 
so much as to produce non-vanishing higher Fourier components, there is a relatively long 
tail to the distribution which may yield non- negligible effect. 

For small smearing, xi,X2 ~ 1, the polar angle dependence of the angular smearing 
can be shown to go as: 

<5^iab^^^Y^sin0. (2.18) 

As expected, the angular smearing is the largest when 6 = tt/2. 

The energy fraction is smeared by about 0.01, which is comparable to the calorimeter 
resolution. This is again accompanied by a long tail. As the cross section is proportional 
to (y^)", we expect that this suppression in the effective centre-of-mass energy may lead 
to a significant suppression of the cross section. 



3. Results 

In this Section, we follow our analysis numerically. The parameters related to the graviton 
emission, the fundamental scale and the number of extra dimensions, are taken to be 
Md = 1 TeV and n = 2 for now. 



- 12 - 



In Fig. the cos 9 dependence for the process e+e~ — > Gj at y/s = 500 GeV is shown, 
together with the corresponding plot for on-shell production. A cut on the energy 
fraction of the photon is taken as x-y > 0.05. In the plot for on-shell production, we also 
show the distribution for the production of a hypothetical scalar Z^. The normalization 
is arbitrary in this case. For this set of parameters, we see that the total rates are quite 
similar. 




Figure 4: The subprocess level cos distribution, for the graviton events integrated over 0.05 < 
Xj < 1 (left), and the on-shell — > i/p events (right). For the events, we also show the 
hypothetical production of a spin-0 object of the same mass (lower curve). The normalization is 
arbitrary in this case. 



From Fig. we see that at first sight, there is no clear distinction between the distri- 
butions due to the production of objects of different spin. We therefore look mainly at the 
moment distribution, as explained in the preceding Sections. 

We then apply the Fourier expansion to this moment distribution, as expressed in 
Eqns. ( p. 7] ) and (2.8). As mentioned before, this is not the only method to extract the 
Fourier coefficients, and it may turn out that a based fit is more practical. We take 
this approach because it is the most theoretically elegant, and also because it provides a 
convenient way of parameterizing the two dimensional distribution of x-y and cos 6. 

Fig. ^ shows the Fourier coefficients for graviton events, before applying any correc- 
tions, as functions of the energy fraction x^. We see that only the first three coefficient 
functions are non-zero, as expected. 

When we truncate the distribution to account for the finite coverage of the calorimeter 
as stated in Sect. 2.4, am for large m acquires an artefactual contribution and becomes 
non-zero. However, we recover the original distribution by the extrapolation procedure 
proposed above. The situation is illustrated in Fig. ^. 

As the corrected distribution shown in Fig. ^ looks almost identical to the original 
distribution shown in Fig. we conclude that the procedure is sufficient for our purpose. 

Let us now take into account the effects of bremsstrahlung and beamstrahlung and, 
perhaps less importantly, the calorimeter energy resolution. For the beamstrahlung and 
calorimeter energy resolution, although the choice is arbitrary and leads to little difference 
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Figure 5: The Fourier coefficients in units of fb GeV^ for graviton events, before applying any 
corrections. 
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Figure 6: The Fourier coefficients for graviton events, after including the calorimeter coverage 
(left), and after the correction by extrapolation (right). 



in the result, let us adopt the TESLA parameters. 

The effect of the ISR is explained and illustrated in Sect. |2.5| . Applying the same ISR 
spectrum to the signal, we obtain the results shown in Fig. |^. We see that the coefficients 
have dropped significantly, and there is further suppression as 1. This is as expected. 

The coefficient as remains zero even with the inclusion of ISR. 

Let us now turn our attention to the SM continuum background. The distribution 
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Figure 7: The Fourier coefRcients for graviton events, including all corrections. 



is shown in Fig. |^. Again, the rate is quite similar to the signal for this set of extra 
dimensional parameters. The dominant part of the background distribution for < 0.8 
is due to the mediated contribution. The rise of the cross section near x-.y = 1 is 
due to the Z^. The coefficients 02 and 03 do not have a peak for this reason. However, 
there is the interference term between the and mediated contributions, and this 
interference term gives rise to the observed behaviour near = 0.95. With the inclusion 
of ISR, the peak broadens. At the same time, the angular smearing due to the ISR gives 
rise to additional contributions to 02 and 03 at the Z^ peak, which lead to the interesting 
behaviour observed near that region. 

From these figures, it is clear that the proper simulation of the ISR is essential to 
understanding both the signal and the background. 



As noted in Sect. 2.3, 02 and higher coefficients are in general non-zero for ^/s suf- 
ficiently higher than Mw. The explicit evaluation, as shown in Fig. ^, confirms this. 
Although these coefficients are non-zero, we note, again, that the shape of the background 
is sufficiently well understood that the distribution can be subtracted away, either from 
the bare distribution or the Fourier coefficients. 

Finally, let us turn our attention to the centre-of-mass energy, y^, dependence of 
the distributions, and the dependence of the signal distribution on the number of extra 
dimensions n. 

In order to understand the centre-of-mass energy dependence, we calculated the signal 
and background Fourier coefficients, which have so far been calculated for ^/s = 500 GeV, 
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Figure 8: The Fourier coefficients for tlie SM continuum background, witlrout the experimental 
and bremsstrahlung corrections (left), and including all corrections (right). 

for ^/s = 1 TeV, including all corrections. The result is shown in Fig. |9[ 

/fb GgV^ /fb GgV^ 




Figure 9: The signal (left) and background (right) Fourier coefficients, at y/s = 1 TeV, including 
all corrections. 



We see that there is a significant change in the shape of the background distribution 
when the centre-of-mass energy is raised. The signal distribution does not change much, 
but if there is any difference, it is expected to be due to the increased beamstrahlung. The 
total rate is increased for the signal and reduced for the background. For the case of the 
background, although the cross section is reduced, the size of the moment distribution 
is increased because of the centre-of-mass energy dependence of the total effective e~e~^ 
charged current cross section which, in the limit of large energy, becomes flat. 

In order to understand the dependence on the number of extra dimensions n, we 
calculated the signal Fourier coefficients for n = 4 and n = 7, at y/s = 500 GeV, again 
including all corrections. The result is shown in Fig. IC. 

We see that as n becomes larger, the Fourier coefficients are more suppressed towards 
larger because of the factor (1— x-y)"/^~^. The higher Fourier coefficients, which rise with 
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Figure 10: The signal Fourier coefficients for n = 4 (left) and n — 7 (right). 

x-y as powers of x^, become suppressed, and the moment distribution becomes more 
flat. In this limit, it becomes impossible to determine the spin from the angular distribution 
only, and one would have to resort to the fitting, for example, of the leading Fourier 
coefficient ao as a function of x^ in order to discriminate the extra dimensional hypothesis 
from other hypotheses. This method would also be useful when there is not sufficient 
statistics to utilize the angular distribution. 

It may be argued that the presence of a quartic polynomial in the coefficient oq, as 
shown in Tab. ^ indicates the spin-2 nature. This is less rigourous, and there would remain 
the possibility that the shape is due to the distribution of gravitons. We also note that the 
size of the quartic term is not large. 

To some extent the difference in the shape of the Fourier coefficients can be exploited to 
'measure' the number of extra dimensions, or more precisely the distribution of the Kaluza- 
Klein modes. One can alternatively measure the ^/s dependence of these coefficients, 
proportional to for n extra dimensions, to yield the same information. 

Having outlined the application of our procedure thus, one may like to question the 
statistical efficacy of our observables. Although a formal discussion would require a more 
detailed study, it is a simple matter to have a rough estimation of the statistics. If we 
assume a relatively low integrated luminosity C = 100 fb~^, and look, for instance, at the 
numbers at 500 GeV shown in Fig. |^, and take 0.8 < x^ < 0.9, say, a quick estimation of 

the numbers show that the significance defined as 'S'/\/^'= \J C^x^a^j (E^ao) ~ 2 so that 
95% confidence level can be attained even with low luminosity for this set of n and Mo- 

4. Conclusions 

We have proposed a way of measuring the 'spin' of invisible objects, including gravitons, 
produced at future LC. 

Our method is to incasiirG tli6 Fourier coefficients Quii^x-y ) (m = 0,1,2- ••) of the E^ 
moment distributions. The presence of terms up to as and the absence of the higher 
coefficients indicates the spin-S nature of the invisible object. 
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We have provided a theoretical justification for this observation, and have carried 
out a simulation based on this principle. Our simulation includes the consideration of 
beamstrahlung, bremsstrahlung, calorimeter coverage, calorimeter resolution and the back- 
ground, and we believe that the simulation is realistic. The ISR, composed of beam- 
strahlung and bremsstrahlung, has a particularly marked effect on both the signal and the 
background. Our study is nevertheless incomplete, as the statistical study and the error 
estimation have not been performed. When doing so, it is desirable that further attention 
is paid to the effect of the hard/multiple bremsstrahlung contribution which must be dealt 
with more carefully. 

This method as stated is useful when the number of extra dimensions n is 2. If we 
believe in the astrophysical constraint, this case may be ruled out. For greater n, the direct 
application of this procedure is not suitable because the higher Fourier coefficients become 
small. Even in that case, the method provides a parameterization of the cross section which 
separates the angular and energy dependence. The resulting x^y dependent coefficient ao, for 
example, can be subjected to analysis in order to test the extra dimensional hypothesis. 

We also note that the shape of the Fourier coefficient functions will provide a measure 
of the 'number' of extra dimensions. More generally, the measurement of spin aside, what 
our method provides is a framework by which we may parameterize single photon events 
at future LC. 
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